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In the weak backscattering limit, point contact tunneling 
between quantum Hall edges is well described by a Poissonian 
process where Laughlin quasiparticles tunnel independently, 
leading to the unambiguous measurement of their fractional 
charges. In the strong backscattering limit, the tunneling is 
well described by a Poissonian process again, but this time 
involving real electrons. In between, interactions create es- 
sential correlations, which we untangle exactly in this Letter. 
Our main result is an exact closed form expression for the 
probability distribution of the charge N{t) that tunnels in 
the time interval t. Formally, this distribution corresponds to 
a sum of independent Poisson processes carrying charge ue, 
2ve, etc., or, after resummation, processes carrying charge e, 
2e, etc. In the course of the proof, we compare the integrable 
and Keldysh approaches, and find, as a result of spectacu- 
lar cancellations between perturbative integrals, the expected 
agreement. 

PACS: 72.10.-d, 73.40. Gk 



The study of dynamic fluctuations is one of the most 
interesting aspects in the physics of mesoscopic conduc- 
tors. Most theoretical and experimental effort so far has 
concentrated on shot-noise [|l|: shot noise is a conse- 
quence of the quantization of charge, and can be used 
to obtain information on a system, that would not be 
available through standard conductance measurements. 
One of the most spectacular progress in this direction 
for instance has been the measurement of the charges of 
Laughlin quasiparticles in fractional quantum Hall de- 
vices 0-0]. In the ideal experiment of interest here, 
one considers tunneling between two quantum Hall edges 
through a point contact. In the weak-backscattering 
limit, a picture where Laughlin quasiparticles tunnel in- 
dependently becomes exact to first order in the current, 
leading to a Poisson statistics for carriers of charge ve. In 
the strong backscattering limit meanwhile, a similar pic- 
ture holds but this time for electrons of charge e [^|-0 . 
The complete noise in between these two regimes has 
been calculated exactly in the field theory limit [§|,||, 
but this does not shed too much light on the statistical 
fluctuations characterizing the tunneling of the charges: 
questions of interest include the nature of corrections to 
the Poissonian limit, and the striking cross-over between 
Laughlin quasiparticles and electrons. 

Study of higher moments of the current is a general way 
to get additional information on the physical processes 



involved in the shot noise: to give trivial examples, a 
vanishing of all higher cumulants would correspond to a 
Gaussian process, while flnding that all the cumulants 
are equal would be indicative of a Poisson process. 

The most natural quantity to consider is the kth cumu- 
lant of the charge N{t) that tunneled between the edges 
during the time t. In terms of the current, this cumulant 
can be expressed as 



(iV,(i)),- / dt^...dt,{I{t,)...I{t,))^ 



(1) 



A knowledge of all the cumulants allows then (under rea- 
sonable regularity assumptions) a reconstruction of the 
probability distribution of the variable N{t) itself. 

So far, cumulants have been studied mostly in the case 
of a pure partition noise, where noninteracting particles 
tunnel with a constant probability r llfl]. In this paper, 
we solve the much more difficult situation of the tunnel- 
ing experiments in P 0, where interactions are crucial. 

The calculation is done at temperature T = 0, and in 
the field theory limit. The first condition corresponds 
essentially to the results reported in ||^-^ ; it is not clear 
however that the scaling regime was reached in these ex- 
periments, and more work in that direction is probably 
needed. Observe that, away from the scaling regime, re- 
sults are nonuniversal, and strongly dependent, for in- 
stance, on the shape of the point contact, hence limiting 
the interest of a theoretical analysis. 

The solution of the tunneling problem in the scaling 
limit was first proposed in [|ll| . It is based on the rather 
well known observation that the quantum field theory 
of interest - the boundary sine-Gordon (BSG) model - 
is integrable [|l2| , but uses a new interpretation of this 
solution in terms of massless scattering |13 and a Boltz- 
mann equation for the integrable quasiparticles. The off- 
shoot is that current and noise can be calculated using a 
model of massless particles with factorized scattering - an 
energy dependent version of the Haldane statistics [Q . 
The current for instance reads / = J^ dO (p+ — p^)t{9) 
where p± are densities of solitons and anti-solitons, 9 the 
rapidity (parametrizing the energy by _E ex e^), t{9) is 
a rapidity-dependent transmission coefficient, whose ex- 
pression is given in B, and we have set e — h — vp = i- 
The characteristic energy scale associated with the im- 
purity is Tb oc e^ . While t follows from the solution 
of the boundary Yang-Baxter equation [n2[ , the densities 



are solution of a set of integral equations known as ther- 
modynamic Bethe ansatz (TEA) equations. The noise 
has a more complicated expression pi, because thermal 
fluctuations of particle numbers at different energies are 
still coupled. However, in the T — limit, only shot 
noise remains. One finds therefore [nl| 



< iVi >=tl = t 



d9p{9)T{9) 



together with 



< TV, >^= t 



dep{e)[r{e)^r'{0)] 



(2) 



(3) 



Here, only solitons contribute, because at T = the 
ground state is made of only one kind of particles, with 
rapidities ranging from — oo (zero energy) to a voltage- 
dependent value of A corresponding to the Fermi energy. 
It was also discovered in |^ that the noise is related 
with the current in a surprising way. 



<r >-- 



di 



2(1 - v) dOi 



(4) 



where what was called < /^ > in [p| corresponds to 
J < N2 >c in the present notations. This relation is rem- 
iniscent of the fluctuation dissipation theorem (which of 
course does not hold in this non-equilibrium situation). 
It occurs technically because the transmission probabil- 
ity T behaves essentially like a Fermi filling fraction, since 
it obeys the basic identity 



mil -Tie)] 
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In the following we introduce the new variable 9b = 

We shall now use the logic of M to compute all the 
other cumulants of N{t) in the tunneling problem. For 
this, suppose first there were only one rapidity, and intro- 
duce for the corresponding quasiparticle a number n = 1 
if it fiips its charge when going through the impurity, 
n = if it does not. One has then, (n) = r. We can 
now consider the higher cumulants of the distribution of 
probability of n. First, recall that 

< 71l >c = < n >= T , 

< "2 >c = < n^ > - < n >^= T -T^ , 

< n3>c = < n^ > -3 < n >< n^ > +2 < n >^ 

^t-3t^ + 2t^ ... . (6) 

From the basic identity (||), we see that < 712 >c= 
IV^^" ■ This result generalizes as follows 



QJ-l 



< rii >c 



(7) 



The proof is quite simple. Consider the Fourier transform 
of the probability distribution of the variable n: p{k) = 
J dnp{n)e^^". Since n takes only two values, it reads 



^ik 



simply p(fc) = l-|-r(e — 1). It follows that dp/dk = ? 
which we can rewrite as dp/dk = iTp + iT{\ — T){e'^^ — 1). 
Use of the basic identity allows us to identify the latter 
term with a derivative with respect to 9'g . Thus we have 



dp 
dik 



dp 

de^ 



p < ni >c 



(8) 
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The cumulants expansion \np = X]i=o "^ ""i ^^ 
then leads us to Eq. (^. 

Now, in the problem of interest, we have charges tun- 
neling at arbitrary rapidities, so the quantities of interest 
are actually integrals. The average of N(t) (the charge 
that went through the impurity in the time interval t) is 
related to the current calculated in [p|. We have 



< N{t) >= t 
and more generally, 
<Nj{t) >c=t 



depi9)T{9)=tI{0B), (9) 



d9p{9) <nj>ci9) . (10) 



If follows therefore that relation (0) also holds for the 
cumulants < Nj >c. The Fourier transform of the prob- 
ability distribution of the variable N obeys then 



-HOb) 



d9B f^, J^- id9By 

^t\l{9B+ik)-I{9B) 



(11) 



Determination of the integration constants leads then to 



lnP(/c) = 2 



1 



d9'B<I 
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Using the expansions for the current derived in |]l6| leads 
to the strong-backscattering expression for the Fourier 
transform of the probability 



P := J]^ exp i yt (e*'=" - 1; 



a„{l/iy) ( V 



n=l 



T' 
^B 



2n{\-v)lv^ 



a^iy) = (-1)"+^ ^^f^^^"'^) , . (12) 

A similar expansion is possible in the weak-backscatter- 
ing limit, giving rise to 



P = e 



vVikt 



(13) 



n°pf"(--"-i)^(i) 



2n(i/-l)" 



In these expressions, T^ is related to Tg (and to the 
bare coupUng constant) in a way made explicit in p6[ . 
Notice that the entire expansions around the two limits 
are completely dual to each other. 

The physical meaning of these expressions is quite fas- 
cinating. To elucidate it, let us recall that for a vari- 
able X taking integer values, this distribution is px — 
m-^ ^j^. Here, m is the average of X. A simple property 
of the Poisson distribution is that all its cumulants are 
equal, so m is the common value of all of them. Equiv- 
alently, P{k) — exp[m (e*'"' — l)]. Suppose now that the 
charge observed in the time interval t were the result of 
a Poisson process for particles of charge one going across 
the barrier, contributing a current /i, plus a Poisson pro- 
cess for particles of charge two contributing a current I2 , 
etc. The final form of the Fourier transform of the prob- 
ability distribution would then be 



Pik) = J] exp [t (e*^-" - 1) /„/n] 



(14) 



This is extremely close to the expression dl3), the only 
difference being that the signs are not quite right: while 
the classical process just considered requires all the /„ to 
be positive, the coefficients in Eq. (n3) are not. Certainly, 
the first contribution is indeed a Poisson process for the 
tunneling of electrons: but the would be tunneling of 
pairs of electrons (and of multiples thereof) comes with 
the wrong sign, a result of quantum interference effects 
we will analyze in more details below. 

The weak-backscattering expansion is of the form 



P{k) 



uVikt 
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l)/„/n . (15) 



The meaning is quite similar: apart from the overall 
factor which corresponds to the current in the absence 
of tunneling (and therefore, absence of fluctuations at 
T = 0), observe that the exponents are now of the form 
g-ifci/n^ The factor v occurs because now we are deal- 
ing with tunneling of Laughlin quasiparticles and mul- 
tiples thereof; the minus sign, because their tunneling 
diminishes the current, instead of building it up as in the 
strong backscattcring limit. The situation for the signs 
is now dependent upon v. Using the inversion identity 
for gamma functions, the sign of /„ for n > 1 is the one 
of cosnTTz/, and therefore positive for the first few values 
of n provided v is small enough. The image of clusters of 
Laughlin quasiparticles tunneling independently with a 
classical Poisson process is therefore quite good for bun- 
dles with modest n, in the small v limit. As v goes to 
zero, all the coefficients are positive, but the fiuctuations 
disappear, as one reaches the classical limit: 



P — exp < vVikt 
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Notice also the remarkable fact that each of the 
"means" of the would be Poissonian distributions are, 
apart from a common factor, pure power laws of the di- 
mensionless coupling constant {V/Tg)'^'-'^~^\ This result 
precludes the /„ from being all positive: in this case in- 
deed, the expansions would have a singularity on the real 
axis, which is nonphysical, except for j/ = 0. 

We stress here that the Poissonian aspect occurs only 
after integration over the rapidity variables, thanks to 
the basic identity (||). At any given rapidity (that is, 
energy), the noise is simply a partition noise, and 



p=l 



T (e*'= 



is not Poissonian, except for very small values of t. Re- 
markably therefore, the interactions have actually made 
the final result simpler than in the noninteracting case! 
The above results rely entirely on the integrable ap- 
proach. We now would like to reconsider them in the 
light of the nonequilibrium Keldysh or Feynman- Vernon 
method. To do so, it is convenient to use the correspon- 
dence of the BSG model with the Schmid model [0 of 
dissipative quantum mechanics. This model describes 
a quantum Brownian particle coupled to an Ohmic heat 
bath and moving in a tilted cosine potential. The strong- 
backscattering limit corresponds to the tight-binding 
(TB) limit of this model. In the TB representation, the 
dimcnsionless viscosity K is related to the filling fraction 
V hy K = l/i/, and the bias energy e/27r corresponds to 
the voltage V. The scaled point contact interaction T^ 
is related to the bare transfer matrix element A and the 
cutoff ojc as T\ 



12-2/u 
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[^^'<lT\K)]^^/u^l'<U 



To proceed, suppose we had a general Poissonian TB 
transport model with transition rates 7^ describing di- 
rect forward (+) and backward (— ) transitions by m 
states. Assuming statistically independent transitions, 
the Fourier transform of the probability distribution 
P{k) = X^n s'''" -P("') ^t t™6 t is found as 
00 

P(k) = n exp[t(e''=" - 1)7+ + t(e-"=" - 1)7,7] ■ (16) 

n=l 

At r = 0, backward moves are absent, and we have 

00 
P{k) = n exp[t(e»^" - 1)7+] . (17) 



The expressions (g^ and (|1^) are equal, but the interpre- 
tation for the terms n > 1 is somewhat different. While 
(nj) describes joint transport of two particles, three par- 
ticles, etc., through a single point contact, the expression 
( |l7| ) represents nonsequential (coherent) tunneling of a 
single particle across one intermediate TB state, two in- 
termediate TB states, etc. We shall now show, using the 
Keldysh or Feynman- Vernon formalism, that this Pois- 
sionian description of the TB model at long times is in 



fact exact. Moreover, it turns out to be possible to cal- 
culate the rates at T = perturbatively, confirming the 
results of the integrable approach at the lowest orders, 
and sheding light on their general structure. 

It is indeed straightforward to derive the exact se- 
ries expression for the generating functional P{k, t) of 
the above dissipative model |lq| . The Laplace transform 
P{k,X) is the analog of an isobaric ensemble of a classical 
gas of charges stringed along the forward and backward 
path. The charge interaction factor for m time-ordered 
charges {uj = ±1} on the forward path q and I time- 
ordered charges {vi = ±1} on the backward path q' is 

^[qm,q'l] = CXp I ^ ^ UjUiQ{tj - u) 
^ i=2 1=1 
I j — 1 I m ^. 

j=2 i=l 1=1 j=l -* 

where Q(t) = 2K\ii[{ujc/2Tr^T)smh2T:^Tt] +iTTKsgnt. 
For a charge sequence contributing to P{n), we have the 
constraints Y], m,- = 2n and V ■ Vi = 2n. 

— ' A 1. — 

The expression ( |16D follows from a cluster decomposi- 
tion of the series for P{k, A — > 0). The clusters are the 
A-independent (irreducible) path sections which start and 
end in diagonal states of the reduced density matrix. By 
definition, the clusters are noninteracting and therefore 
separated by a factor 1/A. Path segments with interme- 
diate visits of diagonal states are reducible, i.e., factorize 
into clusters of lower order. Observe that after subtrac- 
tion of the reducible components always an irreducible 
part is left. The "transition rates" 7^, which do not nec- 
essarily have to be positive, can be identified as the sum 
of all clusters which interpolate between the (arbitrary) 
diagonal state m and the diagonal state mzLn. All the 
charge sequences contributing to the rate can be split 
up in subsets a, 
arrangement {ui, 
Keldysh contour. The property Q(<-z/27rr) = Q*{t) en- 
sures detailed balance for the dynamics, actually already 
for every subset a, 7!"^" = e-^"/2'^^7i"^+ [||. 

At r = 0, the clusters formally appear as the series 



7, 



± _ 



J2a 7" ■ ^"^ ^^'^h subset, the 
j } is fixed and time-ordered along the 
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in which the coefficients /A are given in terms of a sum 
oi2{n + l) — 1-fold integrals, each of them representing a 
particular arrangement of 2(n+Z) charges. Using detailed 
balance for the subsets, 7™ = 0, a multitude of rela- 
tions between the various integrals of the same order can 
be derived |15|. We have analyzed in detail all clusters 
up to order ir . It turns out that 7|j^ is of order x^ , and 
there are no contributions of order a;"* and x^ . Similarly, 
72 is found to be of order a;**, and there is no contribu- 
tion of order x^ . This indicates that there are formidable 



cancellations between the various charge sequences, and 
only those contribute where all charges {uj\ and {vi} are 
positive. From this, it is easy to conjecture, that in the 
series for 7+, at T = 0, only the m — Q term is left. 



7+ = (e/27r)a;2«/("), 



(19) 



in agreement with the result from the integrable ap- 
proach. Namely, a moment relation analogous to (ra), 
< Nj{t) >c=(xd/dxy-^ < Niit) >c, would follow di- 
rectly from ( |l7| ) with (|ig|). A general proof of (pjj ) for 
arbitrary n is still missing: it may well require an analysis 
of the problem similar to the one in Refs. ||l9|, |2C| ]. 

As a final remark, we observe that ( |lj^ ) leads to the re- 
lation fi°^ = [2i-^7r/r(i^)]2"a„(/-i:)/n, and thus allows 
us to construct the full probability distribution for the 
displacement of the quantum Brownian particle too [la . 
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